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The molecular basis of type 2 diabetes predisposition at most 
established susceptibility loci remains poorly understood. KCNQ1 
maps within the lip 15. 5 imprinted domain, a region with an es- 
tablished role in congenital growth phenotypes. Variants intronic 
to KCNQ1 influence diabetes susceptibility when maternally 
inherited. By use of quantitative PCR and pyrosequencing of hu- 
man adult islet and fetal pancreas samples, we investigated the 
imprinting status of regional transcripts and aimed to determine 
whether type 2 diabetes risk alleles influence regional DNA meth- 
ylation and gene expression. The results demonstrate that gene 
expression patterns differ by developmental stage. CDKN1C 
showed monoallelic expression in both adult and fetal tissue, 
whereas PHLDA2, SLC22A18, and SLC22A18AS were biallelically 
expressed in both tissues. Temporal changes in imprinting were 
observed for KCNQ1 and KCNQ10T1, with monoallelic expres- 
sion in fetal tissues and biallelic expression in adult samples. Ge- 
notype at the type 2 diabetes risk variant rs2237895 influenced 
methylation levels of regulatory sequence in fetal pancreas but 
without demonstrable effects on gene expression. We demonstrate 
that CDKN1C, KCNQ1, and KCNQ10T1 are most likely to mediate 
diabetes susceptibility at the KCNQ1 locus and identify temporal 
differences in imprinting status and methylation effects, suggesting 
that diabetes risk effects may be mediated in early development. 
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The translation of established type 2 diabetes risk 
variants into an improved understanding of dis- 
ease pathology is challenging. Progress has been 
made primarily at the few loci where causal alleles 
are coding (1^), but disease mechanisms are unclear for 
the majority of loci mapping outside coding regions. 

The KCNQ1 locus harbors at least two independent 
regions of association with type 2 diabetes risk (intron 10 
and intron 15), both acting through impaired islet function 
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(5-11). KCNQ1 itself encodes the K v 7.1 voltage-gated po- 
tassium channel subunit, which is expressed in human 
p-cells (12) but plays an uncertain role in insulin secre- 
tion. Neither patients with cardiac arrhythmia caused by 
KCNQ1 mutations nor Kcnql -null mice demonstrates 
hyperglycemia or glucose intolerance (13,14), whereas 
KCNQ1 knockdown in human islets does not alter insulin 
secretion (15). 

In accordance with the location of KCNQ1 at the im- 
printed lip 15. 5 region, associated alleles at both signals 
confer disease risk only when maternally inherited (16). It 
has been demonstrated, primarily through studies of the 
syntenic region of mouse chromosome 7, that regional gene 
expression is regulated by differential methylation at the 
promoter of KCNQ1 overlapping transcript 1 (KCNQ10T1), 
a nontranslated antisense RNA that regulates maternal- 
specific expression of downstream genes (17) (Fig. 1). 

Disruption of genomic architecture at the lip 15.5 chro- 
mosomal region has a well-established role in Beckwith- 
Wiedemann syndrome (BWS), a congenital overgrowth 
syndrome often associated with hypoglycemia (18). Fur- 
thermore, unbalanced placental expression of two regional 
genes, PHLDA2 and CDKN1C, is associated with intrauter- 
ine growth retardation (19). We hypothesize that type 2 
diabetes risk may be mediated through disruption of meth- 
ylation and imprinted gene expression within the imprin- 
ted cluster. In this study, we perform the first assessment 
of llpl5.5 regional imprinting in adult human islets and 
fetal pancreas and investigate the effect of risk genotype 
status on DNA methylation and imprinted gene expression. 



RESEARCH DESIGN AND METHODS 

Islet and fetal pancreas isolation and DNA/RNA extraction. Human islets 
of deceased donors of European descent were obtained (with research con- 
sent) from the Oxford Diabetes Research & Wellness Foundation Human Islet 
Isolation Facility (n = 30) and the Human Tissue Laboratory at Lund Uni- 
versity Diabetes Centre (m = 42). Fetal pancreas samples (n = 18) were 
obtained with informed consent and ethical approval from the North West 
Regional Ethics Committee. All islet preparations were >80% pure, with RNA 
integrity numbers >7. Donor and purity details are provided in Supplementary 
Table 1. DNA and RNA were extracted using TRIzol (Life Technologies). 
cDNA synthesis. cDNA was generated through random-primed first-strand 
synthesis from 1 ^.g RNA in accordance with standard protocols, including 
treatment with DNAse I. 

Genotyplng. Type 2 diabetes-associated single-nucleotide polymorphisms 
(SNPs) (rs231362 and rs2237895) were selected as the lead SNPs (strongest 
evidence for association) in each of the independent signals. Reporter coding 
SNPs for imprinting analysis were selected to have the highest possible minor 
allele frequencies, maximizing heterozygous samples capable of differentiating 
mRNA products from homologous chromosomes. Genotyping was performed 
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FIG. 1. Schematic representation of imprinting control [as described in (17)] and type 2 diabetes-associated SNPs at llpl5.5. A: Chromosome 
11:2,450,000-2,960,000. The closed circle at the DMR represents a high level of methylation; the open circle represents a low level of methylation. 
The transcribed sequence is shown in dark gray, and untranscribed is shown in light gray. Arrows indicate direction of transcription. B: Chro- 
mosome 11:2,650,000-2,880,000. The smaller region is distinguished by exonic (boxes) and intronic regions of KCNQ1 (dark gray), a region of 
KCNQIOTI transcription (light gray), and relative positions of top disease-associated SNPs (rs231362, chrll:2,691,471, and rs2237985, 
chrll:2,857,194). All genomic coordinates are b37/hgl9 (graphics not to scale). 



using TaqMan chemistry and SDS2.3 allelic discrimination software (Applied 
Biosystems). All SNPs reached geno typing pass rates of >95%, were present 
in accordance with expected (HapMap Centre d'Etude du Polymorphisme 
Humain) minor allele frequencies, and did not depart from Hardy-Weinberg 
equilibrium. 

Sequencing. Primers were designed using web-based software (http://frodo 
.wi.mit.edu). For cDNA sequencing, primers were designed across exons to 
prevent amplification of residual genomic DNA (gDNA). Samples were am- 
plified using AmpliTaq Gold DNA polymerase enzyme and sequenced on an 
ABI 3700 genetic analyzer machine using the BigDye Terminator vl.l Cycle 
Sequencing Kit (all from Applied Biosystems). Results were analyzed using 
Mutation Surveyor version 3.4 software (SoftGenetics). 

Bisulphite treatment and pyrosequencing. DNA (500 ng) was bisulfite 
treated using the EZ DNA Methylation-Gold Kit (Zymo Research). Treated DNA 
was amplified in duplicate using biotinylated primers and Titanium Taq 
reagents (Clontech) before pyrosequencing on a PSQ 96MA machine (Qiagen). 
Analysis was performed using PyroMark version 2.0 software, including cy- 
tosine bases outside CpG dinucleotides to confirm efficiency of conversion. 
Fragment analysis for indel analysis. PCR was performed using HotStar Taq 
and Q-Solution reagents (Qiagen) with FAM-labeled primers. Products were 
analyzed on an Applied Biosystems 3130x1 machine (Dye set D) and visualized 
using Peak Scanner version 1.0 software (Applied Biosystems). 
Gene expression. Genes were selected for analysis on the basis of inclusion in 
the Geneimprint database (www.geneimprint.com) and expression regulation 
by KCNQlOTVs differentially methylated promoter. Quantitative PCR was 
performed in triplicate on an AB 7900HT machine, using inventoried (where 
available) or custom (designed using Primer Express version 3.0 software) 
TaqMan expression assays and cDNA samples diluted 1:10 in 0.01 mol/L Tris. 
Total expression reactions multiplexed a FAM-labeled test assay and a VIC- 
labeled endogenous control assay (one of RNAseP, HPRT, or PP1 A). Analysis 
was performed using the AAC t method. Allele-specific expression analysis 
compared the amplification cycle number at which VIC and FAM fluorescence 
levels crossed a critical threshold, where VIC/FAM probes annealed over the 
two alleles of a coding SNP. 

Statistical analysis. Methylation and expression in islet samples obtained 
from Lund (n = 42) and Oxford (n = 30) were not significantly different (P s 
0.2, Kolmogorov-Smirnov independent samples test). All islet samples, there- 
fore, were combined for further analysis. Analysis was conducted by linear 
regression of expression or methylation on genotype group, with donor age, 
donor sex, sample purity, and center of origin (for adult islets) or days post- 
conception (for fetal pancreas) included as covariates. Analysis was also 
performed using a Kruskal-Wallis analysis of ranks, which retained signifi- 
cance in all cases. 



RESULTS 

Genes at the llpl5.5 cluster show temporal differences 
in imprinting status. KCNQ1, KCNQIOTI, CDKN1C, 
PHLDA2, SLC22A18, and SLC22A18AS (but not KCNQ1DN 
or TRPM5) were quantiflably expressed in both adult islets 
(n = 72) and fetal pancreas (n = 18). To determine the 
imprinting status for each gene, we used common coding 
SNPs or a coding indel (for CDKN1C only) to distinguish 
between mRNA products from homologous chromosomes 
(Supplementary Table 2). Imprinted or nonimprinted ex- 
pression was indicated by mono- or biallelic expression of 
these SNPs in cDNA from samples that were heterozygous 
in gDNA. 

CDKN1C was monoallelically expressed in both adult 
and fetal tissues (Fig. 2A), as demonstrated in all samples 
heterozygous for the coding indel. Imprinted expression, 
however, was not universal among other genes in the re- 
gion. PHLDA2 and SLC22A18 were expressed biallelically 
in all adult and all fetal samples (Fig. 2B). For SLC22A18AS, 
where there were relatively few heterozygous samples, 
cDNA sequencing demonstrated retained heterozygosity 
at all coding SNPs (Supplementary Fig. 1). KCNQ1 and 
KCNQIOTI showed temporal changes in imprinting status. 
Coding SNPs in both genes showed monoallelic expres- 
sion in all fetal pancreas samples but were biallelically 
expressed in all adult islet samples (Fig. 25). These results 
were confirmed by sequencing cDNA from all heterozy- 
gous samples (Supplementary Fig. 2). 
Type 2 diabetes risk genotype status influences re- 
gional methylation in a developmentally flexible man- 
ner. To test the hypothesis that risk SNPs influence gene 
expression through an effect on DNA methylation, we 
quantified methylation at five representative regions of 
sequence. One assay widely used for the clinical diagnosis 
of BWS (20) maps within the differentially methylated 
region (DMR), and two were located at CTCF and putative 
PLAGL1 transcription factor binding sites within the DMR 



988 DIABETES, VOL. 62, MARCH 2013 



diabetes.diabetesjournals.org 



M.E. TRAVERS AND ASSOCIATES 



Adult Islet gDNA 



Adult Islet cDNA 



1 U UL 









1 





i 1 1 



Fetal 



gDNA 



Fetal pancreas cDNA 



B 



Adult islets 



KCNQ1 



O 

I 
c 
o 
0 

(0 

£ 
a 
x 

0) 

o 
o 

< 



KCNQ10T1 



MO 
ISO 



PHLDA2 



SLC22A18 



l l_ 



l I 



KCNQ1 



KCNQ10T1 
o 

§ 



PHLDA2 



SLC22A18 



0 1 2 0 1 2 

Coding SNP genotype 

FIG. 2. A: Imprinting status of genes in the llpl5.5 cluster CDKN1C. 
Fragment analysis traces demonstrate monoallelic (imprinted) ex- 
pression of CDKN1C in adult islet (top) and fetal pancreas (bottom} 
samples. Left: genomic DNA, with the two size peaks characteristic of 
a heterozygote for dell71APVA highlighted. Right: cDNA from the same 
heterozygous samples. In both cases, only one size peak is evident, 
indicating the presence of mRNA from only one chromosome. Every 
sample heterozygous for dell71APVA at the gDNA level appeared 
homozygous at the cDNA level, indicating monoallelic expression. 
B: Imprinting status of genes in the llp!5.5 cluster KCNQ1 , KCNQ1 OT1 , 



(21,22). Two other assays were designed at the DMR 
boundary on the basis of proximity to numerous CpG 
dinucleotides and disease-associated SNPs (rs231354 and 
rs2283202, respectively). To ascertain risk genotype status, 
DNA from each sample was genotyped for SNPs repre- 
senting the intron 10 (rs231362) and intron 15 (rs2237895) 
signals. 

rs2237895 was nominally associated (P < 0.1) with 
changes in methylation status at three of the sites tested 
and at differing developmental stages (Fig. 3). At the di- 
agnostic DMR region and CTCF binding site, methylation 
was higher (7.3 and 5.6%, P = 0.02 and 0.08, respectively) 
in fetal samples homozygous for the risk allele than in fetal 
nonrisk homozygotes. This effect was not apparent in 
adult islets, where risk genotype had no effect on meth- 
ylation levels at these same sites (P > 0.25). Conversely, 
an effect was seen only in adult tissues at the PLAGL1 
binding site. Here, methylation was 1.6% higher in risk 
genotype homozygotes than in nonrisk homozygotes (P = 
0.006) in adult islets, but no effect was seen in fetal pan- 
creas (P = 0.36). No effects were identified from the top 
disease-associated SNP within intron 10 of KCNQ1 
(rs231362, P > 0.1) (Supplementary Fig. 3) or from either 
SNP at the two candidate assays (P > 0.1) (Supplementary 
Fig. 4). 

Type 2 diabetes risk genotype and total and allele- 
specific gene expression. We identified no relationship 
between rs2237895 risk allele number and total gene ex- 
pression of KCNQ1, KCNQIOTI, CDKN1C, PHLDA2, 
SLC22A18, or SLC22A1AS in the samples available for 
study (P > 0.1 for all genes) (Fig. 4). There was also no 
detectable relationship between rs231362 risk allele num- 
ber and total gene expression (Supplementary Fig. 5) or 
between methylation at any of the five tested sites and 
total gene expression (P > 0.1 in all cases). Power for this 
analysis was reduced by undetermined risk allele parent of 
origin. We also examined the effect of risk genotype on 
allele-specific expression levels but identified no impact on 
the balance of expression between chromosomes for any 
tested gene (P > 0.2). 

DISCUSSION 

There is compelling evidence that diabetes risk at the 
KCNQ1 locus is mediated through a gene with imprinted 
expression (16). By demonstrating that PHLDA2, SLC22A18, 
and SLC22A18AS are biallelically expressed in both adult 
and fetal pancreas and islets, we show that they are un- 
likely to be involved in a proximal molecular mechanism 
for diabetes risk. Likewise, any diabetes susceptibility 



PHLDA2, and SLC22A18. Plots demonstrate flexibility of imprinting at 
the llpl5.5 cluster. The jc-axis represents genotypes of reporter coding 
SNPs within the relevant gene (numbering arbitrary). The left- and 
right-hand clusters contain the two homozygote groups, whereas the 
center cluster contains heterozygous samples. The i/-axis represents 
the ratio of mRNA expression level of a fluorescent probe specific to 
one allele of the coding SNP against the expression level of a differently 
labeled probe specific to the alternative allele. The two homozygous 
groups are therefore expected to have relatively high and low ratios, 
representing substantially more amplification of the probe specific for 
one allele. Under biallelic expression (all adult islet samples and 
PHLDA2 and SLC22A18 fetal pancreas samples), heterozygous samples 
form a central cluster with a ratio of ~1, indicating equal expression 
from homologous chromosomes. Under monoallelic (imprinted) ex- 
pression (KCNQ1 and KCNQIOTI fetal pancreas samples), heterozy- 
gote samples separate into two distinct groups, corresponding to the 
two homozygous clusters and indicating expression from only one 
chromosome. 



diabetes.diabetesjournals.org 



DIABETES, VOL. 62, MARCH 2013 989 



IMPRINTING AT KCNQ1 AND DIABETES RISK 




Diagnostic CpG island 



I I ? 



p = 0.02 



Fetal pancreas 

rs2237895 genotype 

n T2D risk n adult n fetal 
alleles islets pancreas 

20 4 



34 
17 



10 

3 



Adult islets 



Fetal pancreas 



FIG. 3. Methylation according to rs2237895 type 2 diabetes risk genotype. Box plots show the effect of risk allele number (.*-axis) on methylation 
level (i/-axis) separated by tissue type. Boxes represent quartiles; whiskers encompass values within 1.5 times the interquartile range. Risk ge- 
notype is significantly associated with increased methylation levels at the CTCF binding site and Beckwith-Wiedemann diagnostic CpG island in 
fetal samples only and at the PLAGL1 binding site in adult islets only. T2D, type 2 diabetes. 



mechanism working through KCNQ1 or KCNQ10T1 is 
likely to be early in islet development because these 
transcripts are imprinted in fetal pancreas but not in adult 
islets. 

The cyclin-dependent kinase inhibitor CDKN1C (encoding 
p57 mF2 ~), imprinted at both developmental time points, 
emerges as a particularly strong regional candidate. 
CDKN1C is expressed by 30-40% of (3-cells in healthy 
individuals but never concurrently with the Ki67 marker 
of cell proliferation (23). CDKN1C expression is abolished 
in the hyperproliferative pancreatic lesions of focal hyper- 
insulinemia, suggesting a key role in regulating pancreatic 
(3-cell proliferation (23). Loss-of-function mutations in 
CDKN1C cause BWS and hypoglycemia, whereas gain-of- 
function mutations in the PCNA-binding domain have 
recently been shown to cause congenital undergrowth 
(18,24). A higher level of DNA methylation in individuals 
carrying more risk alleles is consistent with a disease 
model of reduced KCNQ10T1 transcription, diminished 
repressive histone modifications (25), and increased 
CDKN1C expression, leading to impaired islet proliferation 
or development. 

We demonstrate that risk genotype status is related to 
DNA methylation in a developmentally variable manner. 
Although the relationship between genotype and methyl- 
ation is statistically stronger in adult than in fetal samples 
(probably attributable to improved power from a larger 
sample size), the magnitude of effect appears larger in fetal 
samples. It is noteworthy that the true allele-specific effect 
size is likely to be underestimated because data were 
obtained by pyrosequencing PCR products amplified from 
both chromosomes. Further work in larger numbers of 
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human islet samples, when they become available, will be 
required to explore in more detail the relationship between 
DNA methylation and expression. 

By use of the largest cohort of human islets currently 
available, we have performed the first assessment of im- 
printing status at lip 15. 5 in human adult islets and fetal 
pancreas. The data provide insights into the complexity of 
imprinting at lip 15. 5, highlighting the necessity of per- 
forming functional studies in relevant tissues and at ap- 
propriate developmental stages. The data have significant 
implications for the molecular mechanism by which as- 
sociated variants in this region exert their effect on di- 
abetes risk. 
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